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ANOTHER PEDAGOGY FOR PURE-INTEGER GOMORY* 
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Abstract. We present pure-integer Gomory cuts in a way so that they 
are derived with respect to a “dual form” pure-integer optimization 
problem and applied on the standard-form primal side as columns, 
using the primal simplex algorithm. The input integer problem is not 
in standard form, and so the cuts are derived a bit differently. In 
this manner, we obtain a finitely-terminating version of pure-integer 
Gomory cuts that employs the primal rather than the dual simplex 
algorithm. 
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Introduction 

We assume some familiarity with integer linear optimization; see [T] for a mod¬ 
ern treatment. We assume that A G Z mxn has full row-rank m, c G Z n , and we 
consider a pure integer-optimization problem of the “dual form” 

z := max y'b 

y'A < c'; (Dz) 

y G Z m . 

The associated continuous relaxation is denoted D. 
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This linear-optimization problem has a non-standard form as a point of depar¬ 
ture, but it is convenient that the dual of the continuous relaxation D has the 
standard “primal form” 


min c'x 

Ax = b; (P) 

x > 0. 


We note that the only linear programs that we directly solve have the form P, 
which is in the appropriate form for treating with the primal simplex algorithm. 

For us, the essence of a pure-integer Gomory (or Chvatal-Gomory) cut is to take 
an inequality a'x < /3 that is valid for the continuous relaxation of a pure-integer 
optimization problem, with a G Z™, and produce the valid cut a'x < |_/3J for 
the pure-integer problem. In the classical presentation of Gomory, the inequality 
a'x < /3 is obtained by rounding down the left-hand side coefficients of a “tableau” 
equation of a standard-form problem, which leads to a valid a'x < f3 because the 
variables are non-negative. In our setting, the integer-constrained variables (on 
our “dual side”) are not non-negative, and our “tableau” equations are on the 
other side (i.e., our “primal side”), so we will have to proceed differently. 

In IJT] we briefly summarize the classical presentation of pure-integer Gomory 
cuts. In 21 we present a new variation of a pure-integer Gomory cutting-plane 
algorithm. Our variation avoids the dual simplex algorithm, and rather precisely 
fits into the well-known column-generation framework based on the primal simplex 
algorithm (commonly used, for example, for Dantzig-Wolfe decomposition and for 
the cutting-stock problem; see [5], for example). Furthermore, it turns out there 
is a certain computational economy which we will explain. In 21 we present an 
example. In (21 we present a finite cutting-plane algorithm based on the cuts of 

m 

An advantage of our set-up and finiteness proof is that it does not rely on the 
lexicographical dual simplex method. In the senior author’s (30+ years) experience 
of of teaching non-doctoral engineering students: (i) already the lexical primal 
simplex algorithm is a topic that many students are challenged to comprehend, 
but eventually learn in the context of proving finiteness for the primal simplex 
algorithm and for establishing the strong-duality theorem of linear optimization 
(i.e., the approach in [5], for example); (ii) a quick explanation of the dual simplex 
algorithm is taken as very technical and somewhat mysterious; (iii) putting these 
two topics together to explain the lexical dual simplex algorithm leaves many 
students behind; so (iv) few students can then absorb the standard (and quite 
technical) finiteness proof for classical Gomory cuts, because it rests on the shaky 
foundation that they have for the lexical dual simplex algorithm. 

We note that throughout (2]and 21) , we deconstruct the derivations and proofs, 
to completely expose the movable parts, rather than seeking to make the presen¬ 
tation as short as possible. 
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1. Classic Gomory 

In the classical presentation of pure-integer Gomory cutting planes (see [3] , for 
example), we seek to solve a standard-form linear-optimization problem P with 
the restriction that all variables are integer. Here we assume that b is an integer 
vector. From a primal basis /3 for the standard-form problem (we use r) for the 
non-basic indices), we derive a Gomory cut from any “source equation” 


x Pi “k ^ ^ Q'i,r]j%rij — x fii (-^i) 

i=i 

having x non-integer. The cut is simply 


x 0i + L a L''jjJ x Vj — L x Pi\ > 
j =i 

which is clearly violated by x = x. Introducing a non-negative slack variable Xk, 
we get an equation 


Xfr + 22 [ a hnA x Vj + x k = i 

j=i 

which, subtracting the equation Ei , can be introduced at the current stage as 

n—m 

'y 1 (Loi,%J — ) Xr lj *k X k = ~ x 0i- 

j= 1 

This new variable Xk is an additional basic variable, but it has a negative value 
x k = — x Pi- Naturally, we proceed to re-optimize by the dual simplex 

algorithm, seeking to regain primal feasibility while maintaining dual feasibility. 
As established by Gomory (see [3]; also see [I] for a presentation more closely 
following the notation used here), this can be realized as a finite algorithm by: 

(i) introducing an integer objective variable xq and associated equation Xq — 
dx = 0; here we note that it is important that c is an integer vector, and 
so Xo is an integer on the set of feasible integer solutions of P; moreover, 
the additional equation implies an additional basic variable, which we take 
as xo and deem it to be the first basic variable: that is, /3 <— (0,/?), now 
an ordered list of to + 1 basic indices from (0,1,2,..., ro}, 

(ii) always choosing a source equation Ei with least i among those with x^ i 
non-integer; here we stress the importance of the objective variable xo 
having index 0 and being the first basic variable at the outset, 

(iii) sequentially numbering added slack variables x n +i, x n + 2 , ■ ■ 

(iv) re-solving each linear-optimization problem after a cut via the lexicographic 
(i.e., epsilon-perturbed) dual simplex algorithm (see [7]). 


4 


Gomory did say in [3]: 

“In these proofs we will use the lexicographical dual simplex method 
described in Section 7. It is not implied that this simplex method 
be used in practice or that it is necessary to the proof. It is simply 
that its use in the proof has reduced the original rather long and 
tedious proofs to relatively simple ones. ” 

This proof has endured in all presentations that we know of (e.g, see m pp. 
215-6]; [2 pp. 165-7]; ® pp. 372-3]; ® pp. 285-7]; [TO] pp. 121-3]; HD PP- 
354-8]), and we do not know Gomory’s “original rather long and tedious proofs”. 
Incidentally, many of these published proofs are lacking a bit in complete rigor, 
including the one of the the second author of the present paper (see m- A clear 
unfortunate aspect of the proof is its delicate set up. 

The classical way of doing Gomory uses the dual simplex algorithm because 
the cut-generation methodology seems wedded to a standard form for the integer 
problem that we wish to solve. In what follows, we derive cuts a bit differently so 
that the dual of the continuous relaxation of the integer problem that we wish to 
solve is in standard form. In this way, we simply add columns to a standard-form 
problem and naturally re-optimize via the primal simplex algorithm. Though just 
a bit more complicated in its derivation than the classical Gomory approach, our 
method can be presented and implemented in a unified and simple manner with 
other column-generation algorithms based on the primal simplex algorithm (in 
particular, Dantzig-Wolfe decomposition, the cutting-stock algorithm, and even a 
presentation of branch-and-bound] see m • Moreover, the actual calculations are 
quite straightforward to carry out (see the example in ]J3]). Finally, we wish to 
point out and emphasize that in carrying out the primal simplex algorithm for P, 
every basis has precisely m elements, even as we add columns. This is in sharp 
contrast to the classical Gomory approach, where each cut increments the number 
of basis elements (as well as appends a slack variable). If classical Gomory were 
to be applied to our formulation Dx (which has m variables in n inequalities), 
putting it into standard form would give us a problem with 2m + n (non-negative) 
variables in n equations. So (dual simplex algorithm) bases would have size n and 
would grow as we add cuts. Because of this, the matrix algebra of each pivot-step 
in our approach is simpler. 

We note that [B] addresses extending our approach to the mixed-integer case. 
Another direction that could be explored is how to lift inequalities to strengthen 
them (see [2], for example). 

2. Gomory another way 

Let us return to approaching the pure integer-optimization problem Dx- In 
what follows, we refer to D (the continuous relaxation of Dx) as the dual and P 
(the dual of D) as the primal. Let /3 be any basis for P. The associated dual basic 
solution (for the continuous relaxation D) is y' := c'pAfj 1 . Suppose that tji is not 
an integer. Our goal is to derive a valid cut for Dx that is violated by y. 
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Let 

b := e l + Apr, 

where r £ Z m , and e* denotes the i-th standard unit vector in R™. Note that by 
construction, b £ Z m . 

Theorem 2.1. y'b is not an integer, and so y'b < [y'b\ cuts off y. 

Proof, y'b = y , (e i + Apr ) = y t + {c'pA^)Apr = 4 r • D 

ez 

At this point, we have an inequality y'b < \_y'b\ which cuts off y, but we have 
not established its validity for Dj. 

Let Hi := Af l e l , the i-th column of Af 1 . Now let 

ui := H.i + r. 

Clearly we can choose r £ Z m so that w > 0; we simply choose r £ Z m so that 

Tk > ~[hki\, for k = 1,... ,m. (1) 

Theorem 2.2. Choosing r £ Z m satisfying m, we have that y'b < \_y'b\ is valid 
for D x . 

Proof. Because w > 0 and y'A < c', we have the validity of 
y'Ap(A^ 1 e l + r) < c'p(A^e l + r), 

even for the continuous relaxation D of Dj. Simplifying this, we have 

y\e l + Apr) < yi + c'pr. 

The left-hand side is clearly y'b, and the right-hand side is 


Vi + c 'pr = yi + y'Apr = y'(e l + Apr) = y'b. 

So we have that y'b < y'b is valid even for D. Finally, observing that b £ Z m and 
y is constrained to be in Z m for Dx, we can round down the right-hand side and 
get the result. □ 

So, given any non-integer basic dual solution y, we have a way to produce a 
valid inequality for Dx that cuts it off. This cut for Dx is used as a column for P: 
the column is b with objective coefficient [y'b\. Taking /3 to be an optimal basis 
for P, the new variable corresponding to this column is the unique variable eligible 
to enter the basis in the context of the primal simplex algorithm applied to P 
the reduced cost is precisely 

y'b - [ V'b\ < 0. 


6 


Observation 2.3. The new column for A is b which is integer. The new objective 
coefficient for c is [y'b\ which is an integer. So the original assumption that 
A and c are integer is maintained, and we can repeat. In this way, we get a 
legitimate cutting-plane framework for Dj - though we emphasize that we do 
our computations as column generation with respect to P. 

There is clearly a lot of flexibility in how r can be chosen. Next, we demonstrate 
that in a very concrete sense, it is always best to choose a minimal r £ Z m satisfying 

HD. 


Theorem 2.4. Let r £ Z m be defined by 


rk = - L hki\, for k = 1,..., m, 


( 2 ) 


and suppose that f £ Z m satisfies r < f. Then the cut determined by r dominates 
the cut determined by r. 

Proof. It is easy to check that our cut can be re-expressed as 

yi < L Vi} + (4 ~ y' A t>) r - 

Noting that — y'Ap > 0 for all y that are feasible for D, we see that the strongest 
inequality is obtained by choosing r £ Z m to be minimal. □ 

3. Example 


In this section, we present an example which illustrates the simplicity of the 
calculations. Throughout, we choose r £ Z m to be minimal, as defined in [2] 

Let 

7 8 —1 1 3 \ / 26 

5 6 -1 2 1 J ’ \ 19 

and c = ( 126 141 -10 5 67 ) . 

So, the integer program Dj which we seek to solve is defined by five inequalities 

in the two variables y\ and j/ 2 - For the basis of P, /3 = (1,2), we have 


A = 


An = 


7 8 
5 6 


and hence A a 1 = 


-4 


-5/2 7/2 

It is easy to check that for this choice of basis, we have 


xp = 


2 

3/2 


and for the non-basis r) = {3,4, 5, 6}, we have c' v = ( 5 1/2 1 ), which are both 

non-negative, and so this basis is optimal for P. The associated dual basic solution 
is 

y 1 = ( 51/2 —21/2 ) , and the objective value is z = 463 1/2. 
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Because both y\ and j /2 are not integer, we can derive a cut for Dx from either. 
Recalling the procedure, for any fraction yi, we start with the i -th column Hi of 
H := Ap 1 , and we get a new A.j := e l + Apr. That is, 



In fact, for this iteration of this example, we get the same cut for either choice of 
i. To calculate the right-hand side of the cut, we have 

3 ) = 70 ^ 


y'b = ( 51/2 -21/2) 


so the cut for Dx is 

42/-j + 3 2/2 < 70. 

Now, we do our simplex-method calculations with respect to P. The new column 
for P is A. e (above) with objective coefficient cq := 70. 

Following the ratio test of the primal simplex algorithm, when index 6 enters 
the basis, index 2 leaves the basis, and so the new basis is /3 = (1,6), with 


Ap = 


7 

5 


4 

3 


with objective value 462, a decrease. At this point, index 5 has a negative reduced 
cost, and index 1 leaves the basis. So we now have /3 = (5, 6), which turns out to 
be optimal for the current P. We have 

y' = ( 131/5 —58/5 ) , and the objective value is 2 = 460 4/5. 


We observe that the objective function has decreased, but unfortunately both yi 
and 2/1 are not integers. So we must continue. We have 

^ = ( 1 3 ) ’ and hence A ^ 1 = ( ^ “/f ) • 

We observe that the objective function has decreased, but because both y-\ and 
2/2 are not integers, we can again derive a cut for Dj from either. We calculate 


H. i 




1 

0 



3 

1 




=: A. 7 



—: A.g . 
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Correspondingly, we have y'A .7 = 96 1/5 and y'A.g = 55 2/5, giving us C 7 := 96 
and c§ := 55. So, we have two possible cuts for Dx: 

5yi + 3?/2 < 96 and 3yi + 2yi < 55. 

Choosing to incorporate both as columns for P, and letting index 8 enter the 
basis, index 5 leaves (according to the primal-simplex ratio test), and it turns out 
that we reach an optimal basis /3 = (8,6) after this single pivot. At this point, we 
have 

y' = ( 25 —10 ) , and the objective value is 2 = 460. 

Not only has the objective decreased, but now all of the fji are integers, so we have 
an optimal solution for Dx- 

We wish to emphasize that to take this example with n = 5 inequalities in 
m = 2 unrestricted variables and put it into standard form, we would end up with 
2m + n = 9 variables and n = 5 equations. So, the initial basis for applying the 
classical Gomory algorithm would have n = 5 elements, and subsequent bases after 
cuts would be even larger. In contrast, our bases have m = 2 elements throughout, 
thus making the matrix algebra less burdensome. 


4. Finite convergence 

To make a finitely-converging algorithm, we amend our set-up a bit: 

(i) we assume that the objective vector b is integer, and we move the objective 
function to the constraints; 

(ii) after this, we lexicographically perturb the resulting objective function. 
So, we arrive at 

max y 0 + y'e { i, m] 

2/o - y'b < 0; 

V'A < c'; (D|) 

2/0 € 2 / 

y e z™, 

where e [M] := (C, e* +1 ,..., e- 7 )', and e is treated as an arbitrarily small positive 
indeterminate — we wish to emphasize that we do not give e a real value, rather 
we incorporate it symbolically. We note that if ( 2 / 0 , 2 /) is optimal for Df, then 
y is a lexically-maximum solution of Dx; that is, y is optimal for Dx, and it is 
lexically maximum (among all optimal solutions) under the total ordering of basic 
dual solutions induced by e? 2/i- 
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The dual of the continuous relaxation of D^ is the rhs-perturbed primal problem 


x 0 


= 1 ; 


- bx 0 + Ax = e [lim] ; 

xo > 0 ; 


x > 0, 


Next, we observe that DJ is a special case of 


(P £ ) 


2 := max y e [0jTrl _ 1] 

y'A < c'; 

y e z m , 


(lex-Dx) 


which has as the dual of its continuous relaxation the rhs-perturbed primal problem 

min c'x 

Ax = e [0 , m _i]; (lex-P) 

x > 0. 

So, in what follows, we focus on lex-Dx and lex-P. 


4.1. First pivot after a new column 


The primal simplex algorithm applied to the non-degenerate lex-P produces a 
sequence of dual solutions y' with decreasing objective value y'e [0tm _ 1] . This can 
be interpreted as a lexically decreasing sequence of y. We wish to emphasize that 
after we add a new column to lex-P, on the next pivot (and of course subsequent 
ones), the basic dual solution y lexically decreases. We want to show more. 


Lemma 4.1. If we derive a column from an i for which yi is fractional (in the 
manner o/ Jl|), append this column to lex-P, and then make a single primal-simplex 
pivot, say with the l-th basic variable leaving the basis, then after the pivot the new 
dual solution is 


y = y 


[Vi\ -Virj 

i —;—■’ 

hu + n 


where Hi. is the l-th row of Ag 1 


Proof. This is basic simplex-algorithm stuff, y is just y plus a multiple A of the 
l-th for of A(f . The reduced cost of the entering variable, which starts at \jji\ — yi 
will become zero (because it becomes basic) after the pivot. So 


(U/d - Vi) - A (hu + n) = 0, 


which implies that 

A = L Vi\ - yi 

hu + n 


□ 
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Corollary 4.2. If we derive a column from an i for which iji is fractional (in the 
manner of §2\), choosing r £ Z m to be minimal (i.e., satisfies append this 

column to lex-P, and then make a single primal-simplex pivot, then after the pivot, 
either (yi ,..., iji-i) is a lexical decrease relative to (iji ,..., jh-i) or iji < |j/ij • 

Proof. A primal pivot implies that we observe the usual ratio test to maintain 
primal feasibility. This amounts to choosing 


l := argmin 

l ■ h H +ri>0 


hu + n 


Also, we have 


<o 


lH = Vi 


-*hu 
hu + n 


>0 

Assume that (iji ,..., yi-\) is not a lexical decrease relative to (iji ,..., iji-i ). Be¬ 
cause y is lexically less than y, we then must have hu > 0. 


- , {.Vi] Vi u - _i_ /1 - i - \ 

Vi = Vi + — - r— h li = Vi + (L Vi] - Vi) 


hu + n 


hii 


hu 


■ n 


<yi + (bd - Vi) = \jji\- 


< o 


> l? 


To finish the proof, we need to justify 


hu 

hu + n 


> 1 . 


($) 


A sufficient condition for $ to hold is r; < 0 and hu >0. Taking r to be minimal, 
we have hu + n = hu — \ hu\ > 0 which, together with hu > 0, implies that hu > 0 
and n = - \hu\ <0 □ 

Observation 4.3. We note that we are using the fact that we choose r to be 
minimal to get $ to hold. However, it is not necessary that we choose r £ TP 1 to 
be minimal for the conclusion of Corollary 14.21 to hold. We simply need to have 
ri < 0 and hu > 0 to ensure that $ holds. 


4.2. A FINITE COLUMN-GENERATION ALGORITHM FOR PURE INTEGER-LINEAR 
OPTIMIZATION 

Next, we specify a finitely-converging algorithm for lex-Dx- We assume that the 
feasible region of the continuous relaxation D of Dx is non-empty and bounded. 
Because of how we reformulate Dx as lex-Dx, we have that the feasible region of 
the associated continuous relaxation lex-D is non-empty and bounded. 
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Algorithm 1: Column-generation for pure integer-linear optimization 

(0) Assume that the feasible region of lex-D is non-empty and bounded. Start 
with the basic feasible optimal solution of lex-P (obtained in any manner). 

(1) Let y be the associated dual basic solution. If yi £ Z for all * € 1, then 
STOP: y solves lex-Dx- 

(2) Otherwise, choose the minimum i £ X for which iji ^ Z. Related to 
this i, construct a new variable (and associated column and objective 
coefficient) for lex-P in the manner of choosing r to be minimal. Solve 
this new version of lex-P, starting from the current (primal feasible) basis, 
employing the primal simplex algorithm. 

(a) If this new version of lex-P is unbounded, then STOP: lex-Dx is 
infeasible. 

(b) Otherwise, GOTO steplU 

Theorem 4.4. Algorithm 1 terminates in a finite number of iterations with either 
an optimal solution of lex-Dx or a proof that lex-Dx is infeasible. 

Proof. It is clear from well-known facts about linear optimization that if the al¬ 
gorithm stops, then the conclusions asserted by the algorithm are correct. So 
our task is to demonstrate that the algorithm terminates in a finite number of 
iterations. 

Consider the full sequence of dual solutions y l (t = 1,2,...) visited during the 
algorithm. We refer to every dual solution after every pivot (of the primal-simplex 
algorithm), over all visits to step I2bl This sequence is lexically decreasing at 
every (primal-simplex) pivot. We claim that after a finite number of iterations of 
Algorithm 1 , y l is an integer vector upon reaching step[TJ whereupon the algorithm 
stops. If not, let j be the least index for which ijj does not become and remain 
constant (and integer) after a finite number of pivots 

Choose an iteration T where y T of step [T| has y\ constant (and integer) for all 
k < j and all subsequent pivots. Consider the infinite (non-increasing) sequence 
<Si := yj, yj +1 ,yj +2 , • • •. By the choice of j, this sequence has an infinite strictly 
decreasing subsequence S- 2 - By the boundedness assumption, this subsequence 
has an infinite strictly decreasing subsequence 1 S 3 of fractional values that are 
between some pair of successive integers. By Corollarv l4.21 between any two visits 
to step [l] with yj fractional, there is at least one integer between these fractional 
values. Therefore, S 3 corresponds to pivots in the same visit to step[2bj But this 
contradicts the fact that the lexicographic primal simplex algorithm converges in 
a finite number of iterations. □ 
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